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Abstract. We consider the magnetic helicity balance for the galactic dynamo in the framework of the local dynamo 
problem, as well as in the no-z model (which includes explicitly the radial distribution of the magnetic fields). When 
calculating the magnetic helicity balance we take into account the redistribution of the small-scale and large-scale 
magnetic fields between the magnetic helicities, as well as magnetic helicity transport and diffusion due to small- 
scale turbulence. We demonstrate that the magnetic helicity flux through the galactic disc boundaries leads to a 
steady-state magnetic field with magnetic energy comparable to the equipartition energy of the turbulent motions 
of the interstellar medium. If such flux is ignored, the steady-state magnetic field is found to be much smaller than 
the equipartition field. The total magnetic helicity flux through the boundaries consists of both an advective flux 
and a diffusive flux. The exact ratio of these contributions seems not to be crucial for determining the strength 
of the steady-state magnetic field and its structure. However at least some diffusive contribution is needed to 
smooth the magnetic helicity profile near to the disc boundaries. The roles of various transport coefficients for 
magnetic helicity are investigated, and the values which lead to magnetic field configurations comparable with 
those observed are determined. 
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• i— i , 1. Introduction Diamond (1995), who suggested that the dynamo growth 

/\, , of the large-scale galactic magnetic field will saturate when 

The magnetic fields of spiral galaxies arc thought to be . n , , , , . n , 

J3 , ill, • i .11 its energy is lower than trie equipartition value by many 

"■^ generated by hydromagnetic dynamo action, based on the , r ., , . ., , , , 

f . . . ....... orders oi magnitudes. As it has been recognized quite 

mint action oi the mean helicity of the interstellar turbu- ,, nz , . , , onnn nl , ^. , , nnn\ 

J ; , „ ., . recently (Kleeorin et al. 2000; Blackman & Field 2000), 

lence and the differential rotation (see e.g., Ruzmaikm et ,, . , , . n . , , 

, . , i their arguments arc based implicitly on local conservation 

al. 1988). The growth of a dynamo generated magnetic t v v i i- ■ i 

; ., . .11 !• °t magnetic helicity: as magnetic helicity is a conserved 

held must eventually be constrained by some nonlinear ., , , . r , , . a r 

, • i r i • i quantity, it can only change it there is a liux ot magnetic 

effect. The most simple and straightforward idea concern- , ,. ., , , , . ,-, , . , . r 

. . . . . • • , helicity through the boundaries, or through resistive cl- 
ing the nonlinear constraint is based on equipartition be- r , , . , , , „ , , , , , n 
° nil lects which are, however, very slow. On the other hand, 
tween the kinetic energy of interstellar turbulence and , n , n . , r ,. , ,. ., r 
°^ dynamo models predict growth ol the magnetic helicity ot 
the magnetic energy. More precisely, the a-effect based , , , . c . , , , ,. , 

,,,, .ii ii the large-scale magnetic held, bmcc the total magnetic he- 

on the hydrodynamic helicity is presumed to be quenched ,. ., . , ln , n n , n , 

. .. i - nil Aii i i hcity is locally conserved, the small-scale and large-scale 

by equipartition strength magnetic holds. Although the , .. . .. , ,. .. , , 

. , . ^ ,n. • ,• iir contributions to the magnetic helicity must have opposite 

idea is very crude, this assumption results m models ol . , r . u . „ „ , , , , . . 

. J .A , , , . , . . ., , , signs. Vamshtem & Cattaneo s conclusions follow, since 
galactic magnetic fields which are quite similar to the ob- n i . • n i i j_- ^ 
° ° , . n . , , the small-scale magnetic held must possess magnetic he- 
served magnetic held distributions m spiral galaxies (see ,. . c , . , , c ., . x tv, . a , i 
, , ° _„„ . . v hcity oi equal magnitude (but ol opposite sign). Ihisheld 
Beck et al. 1996, lor a review). , a , , , , , . , , 

' ...... saturates the a-ertect and hence dynamo action, long be- 

This viewpoint was strongly criticized in the early c . .... r , , , , n n i • . • j 

, b J J fore equipartition ot the large-scale held is attained. 

1990s by Vamshtem & Cattaneo (1992) and Gruzmov & ....... , , , . 

Several possibilities lor overcoming the problem have 

Send offprint requests to: I. Rogachevskii been suggested, and new insight into the physical mecha- 
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nisms of nonlinear saturation supports the applicability of 
conventional galactic magnetic field models (for a review, 
see Kulsrud, 1999). In particular, Moss et al. (1999) stress 
that magnetic helicity conservation arguments do not take 
into account additional magnetic helicity generation aris- 
ing from the buoyancy of the large scale magnetic field, as 
suggested by Parker (1992), and showed that galactic dy- 
namo models based on this mechanism yield magnetic field 
configurations that are quite similar to those observed. 

Moreover, Kleeorin et al. (2000) stress that galactic 
dynamo action is impossible without a turbulent flux of 
magnetic field through the boundaries of the galactic disc 
(cf. Zeldovich et al. 1983) and it is quite reasonable to sug- 
gest that magnetic helicity is also transported through the 
disc boundaries. Kleeorin & Rogachevskii (1999) proposed 
a quantitative model for this flux of magnetic helicity and 
Kleeorin et al. (2000) demonstrated that dynamo action is 
saturated at the equipartition level of the large-scale field, 
provided that the magnetic helicity flux is taken into ac- 
count. 

It remains however unclear to what extent the mag- 
netic helicity flux through the boundaries can contribute 
to the resolution of the magnetic helicity conservation 
problem (see e.g. Brandenburg et al. 2001a). In particular, 
saturation processes involve many other effects in addition 
to any flux of magnetic helicity through boundaries. It is 
necessary to make a careful distinction between the helic- 
ities and the a-effect, and to consider nonlinear quench- 
ing of turbulent magnetic diffusivity, etc. A more specific 
point is that Kleeorin et al. (2000) considered a quite crude 
galactic dynamo model, replacing magnetic helicity losses 
by a simple decay term and considering magnetic helicity 
saturation as a process that occurs much faster than the 
dynamo growth. The aim of the present paper is to remove 
some of these restrictions and, to some extent, to clarify 
these questions. 

We first consider a consistent, nonlinear, one dimen- 
sional model for the local galactic dynamo generation, 
based on arguments of magnetic helicity conservation, in- 
cluding turbulent magnetic helicity transport through the 
boundaries, magnetic helicity diffusion due to small-scale 
turbulence and nonlinear a-quenching, and distinguish be- 
tween helicities and the a-effect. Modern theoretical stud- 
ies provide specific forms for the effects under discussion. 
However, we take into account that the forms used are, to 
an extent, uncertain, and could be model-dependent, and 
so we consider our models also in a more phenomenolog- 
ical light, by varying values of parameters by an order of 
magnitude. We obtain a quantitative agreement between 
our model and the observed distribution of galactic mag- 
netic field. Unsurprisingly perhaps, the 'best' values of the 
model parameters differ from the theoretical predictions. 
Note that we do not here take into account any nonlin- 
ear quenching of the turbulent magnetic diffusion, which 
is the subject of a separate study. We believe that the ef- 
fects of these nonlinearities are less important than those 
modifying the a-effect, and so choose here to isolate the 
latter. A quantitative model for a nonlinear quenching of 



turbulent magnetic diffusion has been recently suggested 
by Rogachevskii & Kleeorin (2001). 

Note that Brandenburg et al. (2001b) also obtained 
equipartition field strengths in their 3D simulations with 
periodic boundary conditions. However, the time scales 
were rather large for high magnetic Reynolds numbers. 
The main reason is that the full treatment of the back 
reaction of the magnetic field onto the turbulence also in- 
volves a turbulent diffusivity quenching. Such a quenching 
is beyond the scope of the present paper, and so we can- 
not expect to reproduce this result of Brandenburg et al. 
(2001b) at the moment by using our approach. However 
there is no obstacle in principle in including a quenching 
of the turbulent diffusion in future studies. 

Our result is not completely straightforward regarding 
estimation of the role of the flux through the boundaries 
(cf. Brandenburg et al. 2001a). We find that all the ef- 
fects taken into account participate in the saturation pro- 
cess at a comparable level. Moreover, our model includes 
two kinds of magnetic helicity flux, i.e. explicit magnetic 
helicity transport as well as diffusion of magnetic helic- 
ity through the boundaries due to small-scale turbulence. 
The action of both effects is quite similar, and it is dif- 
ficult to distinguish between them from modelling of the 
steady-state magnetic field distribution. 

We also perform 2D simulations of the galactic dy- 
namo in the framework of the so-called "no-z" model. 
This involves observable properties only, in that it ignores 
the unobservable vertical magnetic field profiles in galac- 
tic discs, but retains the horizontal field structure. The 
steady-state field distributions obtained are quite similar 
to those known from observations, and from models with 
a naive a-quenching. 

2. The equations for the magnetic and current 
helicities 

The equation for the current helicity (b ■ (Vxb)) can be 
obtained using arguments based on the magnetic helic- 
ity conservation law. We introduce the large-scale vector 
potential A, small-scale vector potential a, and the cor- 
responding representations for the magnetic fields, B and 
b. We then write the total magnetic field as H = B + b, 
and the total vector potential as A. = A + a, so that the 
fields are decomposed into mean and fluctuating parts. 
The equation for the vector potential A. follows from the 
induction equation for the total magnetic field H 

dA/dt = vxH -rjVxH + Vip , (1) 

where v = V + U, and V = (v) is the mean fluid veloc- 
ity field, rj is the magnetic diffusion due to the electrical 
conductivity of the fluid and ip is an arbitrary scalar func- 
tion. Now we multiply the induction equation for the total 
magnetic field H by a and Eq. (|l|) by b, add them and 
average over the ensemble of turbulent fields. This yields 
an equation for the magnetic helicity x m — (a ■ b) in the 
form 



N. Kleeorin et al.: The role of magnetic helicity transport in nonlinear galactic dynamos 



3 



d X m /dt+ V • F = -2{uxb) ■ B - 2r)(b-(Vxb)) , (2) 

where F = (2/3)Vx m + (ax(uxB)) - ij(ax(Vx6)) + 
(ax(uxb)) — (bip) is the flux of magnetic helicity. The 
electromotive force for an isotropic and homogeneous tur- 
bulence is 



(uxb) = aB — t\t (VxJ5) 



(3) 



where t]t is the turbulent magnetic diffusivity, and it is 
assumed that a is the total alpha-effect which at the non- 
linear stage includes both the original hydrodynamic and 
the magnetic contributions. Now we take into account that 
X c = t(6- (Vx6)) is proportional to the magnetic helicity: 
X c — X m /{^ nr lTp) (see, e.g., Kleeorin & Rogachevskii, 
1999), where p is the density and r is the correlation time 
of turbulent velocity field. We would like to stress that \ c 
is a modified current helicity, i.e. the current helicity mul- 
tiplied by the correlation time r of velocity field. Thus, 
the equation for \ c m dimensionless form is given by 

+ % + V • F = A{h/l) 2 [R- 1 B-{VxB)- aB 2 } (4) 



T 



(see Kleeorin & Ruzmaikin, 1982; Kleeorin & 
Rogachevskii, 1999). We adopt here the standard 
dimensionless form of the galactic dynamo equations 
from Ruzmaikin et al. (1988); in particular, length is 
measured in units of the disc thickness h, time in units 
of h 2 jryr and B is measured in units of the equipartition 
energy £? cq = ^Airpu*. \ c an d a are in units of a* 
(the maximum value of the hydrodynamic part of the a 
effect). We define R a = ha*/rj T , T = {l/3){l/h) 2 Rm, 
where / is the maximum scale of the turbulent motions, 
Rm = lu*/r) is the magnetic Reynolds number. Also u* is 
the characteristic turbulent velocity at the scale I, p is the 
gas density, and the turbulent diffusivity rjr = lu t /3. For 
galaxies the term \ c /T is very small and can be dropped. 

The turbulent flux F of magnetic helicity is pro- 
portional to the hydrodynamic part of the a effect 
and the turbulent magnetic diffusivity (see Kleeorin & 
Rogachevskii, 1999; Kleeorin et al. 2000). The turbulent 
flux of magnetic helicity serves as an additional nonlinear 
source in the equation for \ c and it causes a drastic change 
in the dynamics of the large-scale magnetic field. Eq. (Q) 
is the basic equation used below to describe the nonlinear 
evolution of the galactic dynamo. We adopt below a spe- 
cific form of the flux F which includes a transport of the 
magnetic helicity, leading to an advected flux through the 
boundaries, as well as a diffusion flux due to the small- 
scale turbulence. 

3. Basic equations 

The standard mean field dynamo equation takes the well- 
known form 



dB 



V x(V x B + aB - r) T V x B). 



(5) 



The analysis leading to the standard local thin disc dy- 
namo equations for an axisymmetric magnetic field can 
be found in, for example, Ruzmaikin et al. (1988). Using 
cylindrical polar coordinates r, 0, z, we obtain the equa- 
tions for the mean radial field B r = R a b r and toroidal 
field B<f, for the local thin-disc af2-dynamo problem as 



db r 
~dt 

dt 



-(a(B)B (t> Y + b'l 
Db r + B'l . 



(6) 
(7) 



Here b' = db/dz, D = r (dfl / dr) h 2 / t\t is the dynamo 
number, and a(B) is the total a effect which is given by 



a(B) 



(8) 



where a" is the hydrodynamic part of the a effect, and a m 
is the magnetic part of the a effect. These quantities are 
determined by the corresponding hclicitics and quenching 
functions, 4> V {B) and 4> m (B), and B = \B\. In particular, 
a" = X V MB), a m = X c (B)^ m (B) and X " = -(r/3)(« • 
(Vxti)). Thus, 



a(B) = X v MB)+X c (B) ( j )m (B) 



(9) 



where the quenching functions 4> V (B) and (j)m{B) are given 

by 



MB) = (4/7)0 m (5) + (3/7)[l-16B 2 
+ 128S 4 ln(l + (8B 2 )- 1 )] , 

<t> m {B) = (3/8B 2 )(l -arctan(V8B)/\/8B) 



(10) 
(11) 



(see Rogachevskii & Kleeorin, 2000). Thus <f) v {B) — 
1/(4B 2 ) and <t> m {B) = 3/(8B 2 ) for B > l/%/8; and 
<j> v {B) = 1 - (48/5)B 2 and (f> m (B) = 1 - (24/5)B 2 for 
B < 1/V8. 

The function x c (B) entering the magnetic part of the 
a effect is determined by the evolutionary equation 



^ = 4^) [B^K-B^-a^B 2 
+ C^(\ X v (z)\ct> v (B)B 2 )}+K(xr , 



(12) 



where we have written B ■ (VxB) = R a (B^b' r — b r B'^), 
and b' = db/dz, and have introduced a coefficient C for 
the magnetic helicity flux transport term, for future con- 
venience. Unless otherwise stated, C = 1. We have also 
introduced an explicit diffusion of x c , with coefficient n. 

This diffusion plays an important role below; however 
it does not follow directly from the arguments of magnetic 
helicity conservation presented in the previous section. A 
consistent way to introduce it into the equation of mag- 
netic helicity balance requires a more detailed analysis of 
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the interstellar turbulence at the nonlinear stage of the 
galactic dynamo. We consider here that the random flows 
present in the interstellar medium consist of a combina- 
tion of small-scale motions, which are affected by mag- 
netic forces resulting in a steady-state of the dynamo, and 
a microturbulence which is supported by a strong random 
driver (supernovae explosions) which can be considered 
as independent of the galactic magnetic field. The large- 
scale magnetic field is smoothed over both kinds of tur- 
bulent fluctuations, while the small-scale magnetic field is 
smoothed over microturbulent fluctuations only. It is the 
smoothing over the microturbulent fluctuations that gives 
the coefficient k. In particular, the relevant analysis can be 
performed using Eqs. (9)— (11) in Kleeorin & Rogachevskii 
(1999). More pragmatically, it is difficult to imagine that 
a turbulent diffusion acts on the magnetic field, but not 
on the current helicity. 

Equations (^])-(|l^) contain the main nonlinearities 
suggested for the nonlinear a effect. For example, the func- 
tion 4> V (B) describes conventional quenching of the a ef- 
fect. A simple form of such a quenching, <j> v = 1/(1 + 
B 2 ), was introduced long ago (see, e.g., Iroshnikov, 1970; 
Riidiger, 1974; Roberts & Soward, 1975). The splitting 
of the total a effect into the hydrodynamic, a v , and mag- 
netic, a m , parts was first suggested by Frisch et al. (1975). 
The magnetic part a m includes two types of nonlinearity: 
the algebraic quenching described by the function <j) m {B) 
(see, e.g., Field et al., 1999; Rogachevskii and Kleeorin, 
2000) and the dynamic nonlinearity which is determined 
by Eq. ( p^ ) . This equation describes the evolution of mag- 
netic helicity, i.e. its production, dissipation and trans- 
port. The governing equation for magnetic helicity was 
proposed by Kleeorin & Ruzmaikin (1982; see also the 
discussion by Zeldovich et al., 1983) for an isotropic tur- 
bulence, and investigated by Kleeorin et al. (1995) for stel- 
lar dynamos, and self-consistently derived by Kleeorin & 
Rogachevskii (1999) for an arbitrary anisotropic turbu- 
lence. Schmalz & Stix (1991) and Covas et al. (1998) also 
investigated related solar dynamo models. Magnetic helic- 
ity transport through the boundary of a dynamo region is 
reported by Chae (2001) to be observable at the solar sur- 
face. The equation for x c m the form of Eq. ( |l2|) was given 
by Kleeorin et al. (2000), but without the term with coef- 
ficient k. Note that the role of a flux of magnetic helicity 
in the dynamics of the mean magnetic field in accretion 
discs was recently discussed by Vishniac & Cho (2001). 



4. The local thin disc model 

4.1. Asymptotic expansions and an equilibrium solution 

We now present asymptotic expansions for a galactic dy- 
namo model determined by Eqs. (p|)-(|l2}). In a steady- 
state, Eq. (|l^) with k = gives 



Then, for fields of even parity with respect to the disc 
plane, Eqs. @ and (0) give 

B'L' + Da(B)B^ = , 



b r = -D^B'l . 



(14) 
(15) 

For the afi-dynamo B w B^. This assumption is justified 
if \D\ > R a , i.e. [Ru] > 1. Thus Eqs. ©-© yield 

[(B') 2 }' - 2CD (\x v (z)\MB)B 2 )' = . (16) 

For an arbitrary profile Ix^O 2 )! an d negative dynamo num- 
ber D, there is an explicit steady solution of this equa- 
tion with the boundary conditions B$(z = 1) = and 
B'^z = 0) = 0, 

B(z) = (C |#|/2) 1/2 / VWW\dz , (17) 

J\z\ 

where z is measured in the units of h, and we have assumed 
that B ^ 1/V8- For the specific choice of the hydrody- 
namic helicity profile |x t, (- 2: )l = sin 2 (7rz/2) we obtain 



J2C\D\ ( 7TZ 

LI, = y. ^5 oq cos( — 



Br 



7T 



B rn cos , 

cq v 2 



(18) 
(19) 



where we have restored the dimensional factor B cq . The 
boundary conditions for B r are B r (z = 1) = and B' r (z = 
0) = 0. The pitch angle of the magnetic field lines is p = 
-arctan(7r 2 /4|i?J) w 14° for \D\ = 10 and R a = 1. Note, 
however, that our asymptotic analysis in the vicinity of the 
point z = 1 is not self-consistent. The numerical results 
(Sect. 4.2) demonstrate that in this region the role of k is 
important. Note also that for profiles with x"( z — 1) = a 
steady-state solution for B r does not satisfy the boundary 
condition B r (z = 1) = 0. This is the reason why we choose 
the profile with X v (z = 1) ^ 0. 

Now we obtain a steady-state solutions for C = and 
k/0. Thus, Eq. ([TJ) yields: 

B^b' r - B'^br ~ a(B)B 2 + k( X c )" = , (20) 

where k = n/A(h/l) 2 . For fields of even parity with respect 
to the disc plane, Eqs. (|), (Q, (|l|) and © give for 
B > l/VS: 

[(B 1 ) 2 2\D\k( X c yy = , (21) 



where B w B^ for the ail-dynamo. Combining Eqs. (|20| 
and (Ell) we obtain 



4(BBW + B'B'") - ('S/Ak)(B') 2 = -\D\[(x v )' 

+ (3/2)( X c )Uo] , (22) 

where B^ IV ^ is the fourth-order ^-derivative. For the spe- 
cific choice of the hydrodynamic helicity profile x v — 
sin(7rz) there is a steady-state solution 

1/2 



a(B)B 2 = B^b' r - B'^br + C[\ X v (z)\MB)B 2 



(13) 



B„ 



B r 



D 



•3/ (8k) 



B eq cos 



1/2 



(23) 



|^|(tt 2 + 3/(8k)) 



£ eq cos — , (24) 
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which exists for the boundary condition (x c )' z=0 = 
— (4tt 2 k + 3/2V" 1 . Note that the term proportional to 
k(x c )" in Eq. ( |l2| ) can be interpreted as a flux of magnetic 
helicity which causes a steady-state solution for the mean 
magnetic field. The pitch angle for this solution turns out 
to be the same as for that with k = 0. 



4.2. Numerical solutions 

We found solutions of Eqs. (|| - |l^) by step by step inte- 
gration, from arbitrarily chosen initial conditions. 

In agreement with our previous results (Kleeorin et al. 
2000), the flux terms connected with nonvanishing coef- 
ficients C and k lead to (large-scale) magnetic field evo- 
lution at significant steady-state magnetic field strengths, 

1. e. at strengths comparable to the equipartition value. 
However, the relative roles of C (magnetic helicity trans- 
port) and k (magnetic helicity diffusion due to the small- 
scale turbulence) appear to be quite unexpected (note that 
the analysis of Kleeorin et al., 2000 was too crude to make 
a distinction between these effects). We obtained steady- 
state solutions even for vanishing C and nonvanishing «, 
while a weak diffusion of magnetic helicity is needed in 
order to get a steady-state with nonvanishing C. This 
fact should be connected with the peculiar behavior of 
the asymptotic solution near to the disc boundaries when 
k = 0, mentioned in Sect. 4.1. However this distinction be- 
tween the role of magnetic helicity transport and diffusion 
seems to be more of a mathematical rather than a phys- 
ical nature. For n = 0, Eq. jl^ ) is a first order equation 
with respect to the spatial variables, whereas it is a sec- 
ond order equation for k / 0, and this difference causes 
problems with the boundary conditions. Presumably, it 
means that the evaluation of the term proportional to C 
in this equation, whilst self-consistent inside the galactic 
disc, needs some improvement near to the sharp bound- 
ary. In general, we find that both effects destroy the local 
conservation of magnetic helicity and that their variations 
affect the final magnetic field distribution in a similar way. 

We used spatial profiles x" = sin 2 (7rz/2), \ v = 
sin(7rz), \ v = sin 2 (7rz), for values of the ratio h/l = 

2, 5, 10, and put C = 1, 10. The solutions are steady, ex- 
cept for D = —1000 and h/l > 5, in the cases with either 
X v = sin(7rz) or \ v = sin 2 (7rz/2), and C = 10. 

Various properties of these solutions are illustrated in 
Figs. 1 - 5. It can be seen clearly that the field strength 
(~ \B<f,\) is typically of order 1 (equipartition) or larger, 
and increases with \D\. The role of the flux transport term 
in Eq. ( |l2] ) can be seen in Fig. 3, i.e., the magnitude of 
the mean magnetic field \B\ increases with the coefficient 
C. 

However, it is clear that a detailed agreement with 
observed galactic magnetic fields needs some effort. In 
many cases we obtain a steady-state magnetic field about 
0.1 — 0.5B cq while observations suggest B w B cq . To get 
this we need \D\ > 20 — 100 (Fig. 2), which is a quite 
restrictive condition, or a value of C that exceeds unity 
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Fig. 1. The z-dependence of solutions for the one- 
dimensional model; dynamo number D = —20, C = 1, 
k = 1 and |x"( 2; )l — sin 2 (7T2:/2). The various curves indi- 
cate results with h/l = 2 (short-dashed), h/l = 5 (dashed) 
and h/l — 10 (solid), a) the toroidal magnetic field B^(z), 
b) the poloidal magnetic field B r (z), c) the total a-effect, 
ot{z), d) the function x c ( z )- 



(Fig. 3). However straightforward theoretical estimations 
suggest C < 1. Even so, the values of C required do not 
appear unreasonable, given the inherent uncertainties of 
the theory. 



5. The no-z model 
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Fig. 2. As Fig. [l], with h/l — 5, and dynamo numbers D = 
-20 (short-dashed), D = -100 (dashed) and D = -1000 
(solid). 

5.1. The model equations 

The basic no-z dynamo model for disc galaxies is described 
in Moss (1995). Phillips (2001) suggested a plausible tun- 
ing of the model and his major amendment, the multi- 
plication by factors 7r 2 /4 of the terms representing the 
z-diffusion of B r and B^,, was implemented here. The no- 
z model differs from the local model of Sect. || in that 
it describes magnetic fields over the entire radial range, 
< r < R, but all explicit dependence on the vertical 
coordinate z has been removed, with z-derivatives being 
replaced by inverse powers of h. The field components B ri 
B<p appearing in the no-z equations can cither be thought 
of as representing mid-plane values, or as some sort of ver- 
tical average of values through the disc. Here we develop 
a version of Eq. (|l2|) for the no-z model, starting from the 
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Fig. 3. As Fig. [j], with D = -20, showing the effects of 
changing C: C = 1 (dashed) and C = 10 (solid); h/l = 5. 



appropriate analogue of Eq. (n2h which is 



R~ X B ■ (V X B) - olB 2 



+A X C - div[V x c ] 
a = x v ®v{B/Jp)+x c ® m {B/yfp) , 



(25) 
(26) 



where V = f2 r is the large-scale velocity (differential 
rotation). The factor p^ 1 multiplying by the right hand 
side of Eq. (|2^) arises because the link between the mag- 
netic helicity \ m and the current helicity x c contains the 
density p, i.e. x c = X m /^ nr lTP- The quenching functions 
<f> v and $ m contain in their arguments the factor 1/yfp 
because they are based on local equipartition, while Eqs. 
( p5| ) and (|2^) presume that the field B is measured in units 
of equipartition at the point r = 0. Note that our model 
includes a flux of magnetic helicity which is not directly 
connected with a mean magnetic energy flux; however an 
energy flux associated with the small-scale magnetic field 
is, of course, possible. 



dt 



The coordinate version of Eq. (25) is given by 

" = P_1 { 4 (^) ^ 1 ^ B 'r- B - B ''p)- aB2 

+(\X V \ $v B 2 )'] + —-[r X v A; 1 B 2 ] + (x c T 



A 2 d_ 

r dr 



dr 



^S-A 2 ^l, (27) 
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Fig. 4. As Fig. |, with D = -20, C = 1, h/l = 5, showing 
the effects of changing the form of x"'- \x v ( z )\ — sin 2 (7rz/2) 
(dashed), and x"( z ) — sin("7rz) (solid). 

where (x c )' = dx c /dz, A = h/R is the aspect ratio, and 
Ap 1 = — |V|p//J is the density stratification scale. 

For the no-z model in the axisymmetric case Eq. ( p7[ ) 
becomes 



at 



[R- 1 (B r B^ + C X v $ v (B/^/p)B 2 



A\2 a 

-«B 2 ] - «X C + — Q- r [r A; 1 B 2 
A 2 9 



H r 

r 9r \ cV 



(28) 



Note that with the no-z formalism, the term B ■ (V x 
_B) vanishes identically. Examination of the solutions of 




0.2 0.4 0.6 0.8 1 




Sect. 4.2 suggests that setting this term to zero may not al- 
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Fig. 5. As Fig. g, withZ? = -20,(7= l,h/l = 5, \x°(z)\ = 
sin 2 (7rz/2), showing the effects of changing k: k — 0.1 
(dashed), k = 1 (short-dashed), k = 10 (solid). The values 
of x c vary in the ranges (-0.5,0.25), (-0.25,0.10) and 
(-0.04,0.15) for k = 0.1, 1 and 10 respectively. 

ways be a good approximation, but experimentation with 
the inclusion of an order of magnitude estimate for the 
omitted term suggests that the no-z solutions are insensi- 
tive to it. We have introduced a coefficient k in the term 
— kx c in Eq. (p8|). The motivation is that this term rep- 
resents the term d 2 x c /dz 2 in Eq. (|27|). The literal ap- 
plication of the no- z ru les would imply k = 1. However 
the results of Sect, [f .2| suggest that in the local model 
X c ~ sin(27rz), giving k w 47r 2 . This is analogous to the 
modification of the representation of the z-diffusion terms 
for the magnetic field suggested by Phillips (2001), and 
mentioned above. Thus we introduce free param- 

eter. The term <$> v (B j J~p) B 2 represents a magnetic 
helicity flux. As it is introduced in a somewhat ad hoc 
manner, we have introduced the coefficient C to represent 
this uncertainty. To determine the magnetic field distri- 
bution along the radius we use a Brandt rotation law, 
fl(r) = fi /[l + (r/^) 2 ] 1/2 with r w = 0.2, and the radial 
density profile p(r) — exp[— (r/r p ) 2 ] with r p = 0.5, so that 
A p = 2r/r 2 . We also set x v (r) = 1. 

One technical point should be noted here, which is im- 
portant when comparing results from the local thin-disc 
model, studied in the earlier parts of this paper, with those 
from the no-z model. For the local thin-disc model, R u = 
rh 2 r]^, 1 dn/dr < 0. By the nature of the model, rdfl/dr is 
the value at a chosen radius in the disc, and r does not 
further occur explicitly in the analysis. However the no- 
z model is global with respect to radius, and the value of 
rdfl/dr varies through the disc, from zero at r — to some 
maximum absolute value; for the Brandt rotation law this 
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value is O.31f2o at r = a/2/ 3r w . For the no-z model the 
global definition is R u = i? w (no— z) — Vl h 2 /r]T > 0. (Less 
importantly, there are also small differences, of order 25%, 
in the effective values of R a occurring in the two approxi- 
mations, even though the formal definitions are the same 
- see Phillips (2001).) To emphasize these distinctions, we 
write the no-z dynamo number as T> = R u (no — z)R a . 

For these reasons, it cannot be expected that there 
will be a precise correspondence between the values of 
the marginal dynamo numbers for these two approxima- 
tions; empirically we found X> cr it ~ 20. We note that 
I -Dcrit | I'D crit ~ 0.4, which can be compared with the maxi- 
mum value of \rdQ / dr\ / = 0.31 for the no-z model with 
our Brandt rotation law. These issues are also touched on 
in Sects. 3.2 and 4 of Moss et al. (2001). 

What seem to be important are the terms with A 2 . 
They represent a diffusion and an advection of helicity 
along the radius, i.e. from the galactic center to its pe- 
riphery. This could lead to a more extended magnetic field 
distribution than given by conventional models. This is 
why the result can be compared with that from the con- 
ventional model with a ~ 1/(1 + B 2 /p). 

5.2. Numerical results 

We integrated the standard axisymmetric no-z equations 
(modified as described above) for B r and B,p, together 
with Eq. (psl), from arbitrarily chosen initial conditions 
for a range of values of C and k, finding steady solutions. 
We present in Fig. 6, as examples which are close enough 
to the available phenomenology of galactic magnetic fields, 
solutions with k = 1, £ = 0, V = 50, C = 3 (dashed line) 
and C — 10 (solid line). 

Again, we can observe that typically \B\ is of order 
unity. The value r — 1 can be taken roughly to corre- 
spond to 15 kpc in dimensional units, so the radial mag- 
netic field reaches its maximal value at a galactocentric 
radius of 3 ... 4 kpc, quite typical for the simplest approx- 
imations of galactic rotation curves. Pitch angles in this 
radial range are about —15° ... — 20°, in agreement with 
available observational data. The values of \ c are m the 
range 1.5 .. . 6, so the back reaction of magnetic field ba- 
sically modifies the dynamo action. The main dependen- 
cies of the results on the model parameters are given in 
Table 1. 

We conclude from the the results shown in Table 1 that 
in the no-z model the helicity diffusivity k does not play an 
important role (in contrast to the ID model of the previous 
section). The only anomalous thing about solutions with 
k = is an extremely long (about 50 diffusion times) 
transient phase when C = 0.1. However this was not a 
slow growth, but rather a growth over a few time units 
followed by a rapid decay and then a slow relaxation. This 
feature was not seen when C = 1, 10, when the magnetic 
fields grew to approximately their final strengths over a 
few time units. 
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Fig. 6. Details of the no-z model, with dynamo number 
V = 50, h/l = 5, A = 0.05, \x v {z)\ = sin 2 (7rz/2). a) the 
toroidal magnetic field -£></, (r), b) the radial magnetic field 
B r (r), c) the pitch angle p, d) the function x c { z )- Broken 
curves are solutions with C — 3, solid with C = 10. 

The interaction of the two flux transport coefficients is 
not straightforward. For C = 0.1, field strengths increase 
with k, whereas for C = 10, they decrease. When C = 1.0, 
the solutions are remarkably insensitive to the value of 
K, independent of the value of T>. These comments apply 
both to solutions with a modestly supercritical dynamo 
number (T> = 50), and a very supercritical V = 400. The 
pitch angle depends predominantly on V - see Table 1. 

6. Discussion and conclusions 

Conventional galactic dynamo models are based usually 
on a naive concept of a-quenching. From the theoretical 
point of view this concept must sooner or later be replaced 
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Table 1. Summary of selected calculations for the no-z 
model. Er/Ecf, is the ratio of energies in the radial and 
azimuthal magnetic fields, (p) and (x c ) are the averages 
with respect to radius of pitch angle (in degrees) and x c 
in the range 0.20 < r < 0.80, and m is the maximum 
value of B^. 



Table 1. 



Ra — 


2, R^ = 


25, V = 


= 50 




C = 0.1 


K 




<P> 


<x c > 







0.089 


-16 


-0.67 


0.028 


0.1 


0.118 


-18 


-0.56 


0.060 


1.0 


0.113 


-18 


-0.42 


0.149 


10.0 


0.103 


-18 


-0.18 


0.320 


C = 1.0 


0, 0.1, 1.0, 10.0 


0.097 


-18 


« 


0.484 


C= 10.0 





0.104 


-18 


6.22 


2.19 


0.1 


0.103 


-18 


6.1 


2.19 


1.0 


0.100 


-18 


5.5 


2.15 


10.0 


0.096 


-18 


2.9 


0.82 




C =- 


1.0 






1.0 


0.13 


-18 


-0.52 


0.073 


10.0 


0.11 


-18 


-0.30 


0.021 


Ra = 


5, R^ = 80, V = 


400 




C = 1.0 


0.1, 1.0, 10.0 


0.01 


-6 


?a 


1.93 



by a better elaborated model for the processes of magnetic 
helicity generation and transport, which lead to dynamo 
saturation. In this paper we have presented a model of 
such saturation, based as far as possible physically on first 
principles rather than on ad hoc parameterizations of tur- 
bulence. This model is shown to reproduce successfully 
basic features of galactic magnetic fields and, when used 
in conjunction with standard galactic dynamo models, can 
be used to obtain detailed models of magnetic field gen- 
eration in specific galaxies. Perhaps not unexpectedly, we 
have to adjust some coefficients of the model in order to 
obtain this agreement. However all terms occurring in the 
model can be justified from first principles. 

Magnetic field distributions obtained with our model 
are compared with that from a model with a conventional 
a-quenching and corresponding values of the control pa- 
rameters, in Fig. 7. It can be seen that results for these 
models closely resemble each other, and for the choice 
C = 3 they practically coincide; certainly their deviation is 
much less than the available accuracy of the observations. 
Note that our model, in the case shown in Fig. 7, is intrin- 
sically quite distant from an a-quenched model. However, 
for C — 1.0 and T> = —50 we find (x c ) ~ 0, so here the 
dynamo saturation must be connected mainly with a con- 
ventional cv-qucnching and n has no influence. We take 
here k — 1 although, as noted, solutions are insensitive 
to n for these parameters. Of course, the last result does 




0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 




-40°- 



0.2 0.4 0.6 0.8 1 

r 



Fig. 7. Comparison of the equilibrium r-profiles for the 
no-z model obtained using two kinds of nonlincarities: the 
nonlinearity considered here for C — 1 (dashed) and C = 
3 (solid), and a standard alpha-quenching a tx 1/[1 + 
B 2 (r)/p(r)] with p = 1 (short-dashed). Notation is as in 
Fig. 6. 



not coincide with that for conventional a-quenching, be- 
cause Eq. ( p8| ) is more complicated than the conventional 
a-quenching, even when (x c ) ~ 0. 

From a practical viewpoint, these results lead to a two- 
fold conclusion. On one hand the research undertaken can 
be considered as a justification of conventional galactic dy- 
namo models based on a naive a-quenching. However the 
physics underlying our model is much more complicated. 
On the other hand it means that the galactic dynamo, 
at least in its standard manifestation considered here, is 
very robust and the resulting magnetic field is practically 
independent of small details of the dynamo model. It fol- 
lows that it is necessary to consider more fragile dynamo 
regimes in order to isolate features specifically connected 
with the details of dynamo saturation mechanisms. 

Our analysis has also shown that the detailed form 
of the flux of the magnetic helicity is not crucial for the 
overall state of the saturated mean magnetic field. The 
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most important factor is to have a nonzero flux of the 
magnetic helicity (e.g., in the form of the transport term 
or in the form of the diffusive term) . 

We note that our results clarify the role of magnetic 
helicity transport in galactic dynamos in the mean-field 
approach only. However, we do not address the question 
of the link between the transport of the total magnetic 
helicity A. ■ H and the total magnetic field H . A recent 
paper by Brandenburg et al. (2001b) on the topic of stellar 
dynamos using direct numerical simulations isolates many 
further questions which are still far from being clear. 

Finally, we note that theoretical ideas do not even ex- 
clude negative values of C. In this case larger steady state 
magnetic fields might be predicted, because the flux trans- 
port could then amplify the mean field and the saturation 
might be connected with the diffusive flux only. However 
the results of simulations with negative C given in Table 1 
demonstrate that the detailed balance between the vari- 
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ous terms in Eq. (28) is too complicated for this simple 
explanation to be valid. 
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